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ABSTRACT 



Molecular FeH provides a large number of sharp and isolated absorption lines that can be used to measure radial velocity, rotation, or 
magnetic field strength with high accuracy. Our aim is to provide an FeH atlas for M-type stars in the spectral region from 986 nm 
to 1077 nm (Wing-Ford band). To identify these lines in CRIRES spectra of the magnetically inactive, slowly rotating, M5.5 dwarf 
GJ1002, we calculated model spectra for the selected spectral region with theoretical FeH line data. In general this line list agrees 
with the observed data, but several individual lines differ significantly in position or in line strength. After identification of as many 
as possible FeH lines, we correct the line data for position and line strength to provide an accurate atlas of FeH absorption lines for 
use in high precision spectroscopy of low mass stars. For all lines, we use a Voigt function to obtain their positions and equivalent 
widths. Identification with theoretical lines is done by hand. For confirmation of the identified lines, we use statistical methods, cross- 
correlation techniques, and line intensities. Eventually, we were able to identify FeH lines from the (0,0), (1,0), (1, 1), (2, 1), (2,2), 
(3, 2), and (4, 3) vibrational bands in the observed spectra and correct the positions of the lines if necessary. The deviations between 
theoretical and observed positions follow a normal distribution approximately around zero. In order to empirically correct the line 
strength, we determined T eff , instrumental broadening (rotational broadening) and a van der Waals enhancement factor for FeH lines 
in GJ1002. We also give scaling factors for the Einstein A values to correct the line strengths. With the identified lines, we derived 
rotational temperatures from line intensities for GJ1002. We conclude that FeH lines can be used for a wide variety of applications 
in astrophysics. With the identified lines it will be possible for example to characterize magnetically sensitive or very temperature 
sensitive lines, which can be used to investigate M-type stars. 
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1. Introduction IWing & Foidl (I1969I) first discovered the molecular band 

around 991 nm in the spectra of the cool dwarf Wolf 359. This 
^. . High resolution spectroscopy of atomic or molecular lines is band was also found fa s . type stars Jy^'j!^ and was iden- 
r- used to measure rotation, magnetic fields, and radial velocity tifie d as the (00) vibrational band of the FeH molecule by 
variations. In solar-like stars, atomic lines are useful to measure Nordh et al ( 1977) An extensive analysis was carried out by 
these quantities since they are numerous, well isolated and suf- p hmips et al ( 1987) They identified seven vibrational bands of 
ficiently narrow. In cooler stars, like M dwarfs, atomic lines are the 4 A _4 A electronic transition of the FeH molecule and pro- 
no longer useful because most of them become very weak, others vided tables wkh molecular constan ts and quantum numbers. An 
become strongly pressure broadened and they are usually over- important theoretical work, partly based on the previous one, 
lapped by strong molecular bands. In these stars molecular lines was carried out by | Du ii cketa iJ (20 03). They computed a line 
are a valuable tool to measure the quantities mentioned above. Hst for the F 4 A _ X 4 A electr0 nic transition and provided large 
However, molecular lines tend to cluster in dense bands e.g. for tables of moleC ular data with quantum numbers and line intensi- 
TiO and VO in the visual spectral range. Only a few molecules t - 

provide absorption lines that are isolated and can be used for _ TX , .,, ,, , . 

\ . ., ■ . • , • FeH absorption bands were also detected in the J- and 

detailed spectroscopic analysis. , T , , • , ,. , . nz — t-. ni ---j 

m_ J; T1 , , ., ^ i i i i_ c //-band with medium resolution spectra (Cushing et al. 2003). 

The FeH molecule provides a particularly large number of T , ... , , „ 4 _ . 4l -. , . .. . . ., , 

j 11 ■ i . j i- • *i. u j / r,c\n mon \ t. In the //-band the E*I 1 - A^Yl electronic transition is visible 

strong and well isolated lines in the z-band (~ 990-1080nm). It rrr . i l Urvink t,, . t- n u j. j . r 

. , ° . ., ... .j. , . , c (Hargreaves et al. 2010). That FeH can be used to determine er- 

is the main opacity contributor in this region for late-type dwarf : — t° — . „ , , i „ . . — n 

, r , , j? ,. , t-tt tective temperatures was shown for example by Schiavon et al. 

stars, and can be used for high precision spectroscopy. FeH pro- rrr^k ,,, — ^ — -nrrri , . r J~T. I 

■j ui j jv *t. * a: ■ .i . 1997) or Wende et al. 2009), and its potential to measure 

vides numerous unblended lines that are sufficiently narrow to . J - ~ — * ~f r , , nr~- m — n 

... , . „ . . / .. . magnetic field strengths was demonstrated by Reiners & Basn 

measure small broadening effects or variations in the line posi- ,~nn^ ™rv-^ t., ■ , , , T'- — — Z r TT~ r 

° (2006, 2007). Theoreti cal work on the magnetic sensitivity of 



X 
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FeH was published bv lAfram et al.1 (1200 7. 2008). 



* Data were taken at ESO Telescopes under the program 79.D- In this paper we use high r esolution spectra o f GJ1002 to 

0357(A) empirically verify the line list of lDulick et al.l (l2003h and identify 
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FeH on a line-by-line basis in the region 989.8 - 1076.6 nm. For 
this purpose we will use Voigt functions to determine the empir- 
ical positions and equivalent widths of the observed FeH lines. 
We identify the observed lines with theoretical ones by hand, 
and confirm this identification using statistical means, cross- 
correlation techniques, and the line strength of the identified 
FeH lines. Furthermore, we correct theoretical Einstein A val- 
ues to account for mismatches in line depth. We also derive ro- 
tational temperatures from the identified lines, and investigate 
under which circumstances they are close to the effective tem- 
perature of the star. 



2. Data 

2.1. CRIRES spectra of GJ1002 

The observational data are CRIRES spectra of the inactive M 
5.5 dwarf GJ1002 (see Fig. Q]). The M dwarf has an assumed 
effective temperature of 3150 K (from the spectral type), and it 
is a very slow rotator (v sin i < 3 km s _I , Reiners & Basri 2007). 
There is also very low H a and X-ray activity from which we can 
assume that the magnetic field strength is relatively low. Due to 
the weak magnetic field and slow rotation, GJ1002 is an ideal 
target for identification of molecular FeH lines. The lines are 
only slightly broadened by the different possible mechanisms in 
contrast to observations in sun spots where FeH lines are always 
influenced by a strong magnetic field. 

CRIRES observations of GJ1002 were conducted in ser- 
vice mode during several nights in July 2007. The entrance slit 
width was set to 0.2", hence, the nominal resolving power was 
R ~ 100000. Four frames with an integration time of five min- 
utes were taken in an ABBA nod pattern for each of the nine 
wavelength settings covering the region between 986 - 1077 nm, 
leaving only one larger gap at 991.15 - 992.45 nm and two 
smaller gaps at 997.15- 997.50 nm and 1057.15 - 1057.65 nm. 

Data reduction made use of the ESOREX pipeline for 
CRIRES. Science frames and flatfield frames were corrected for 
non-linearity and ID spectra were extracted from the individual 
flatfielded and sky subtracted frames using an optimum extrac- 
tion algorithm. The wavelength solution is based at first order 
on the Th-Ar calibration frames provided by ESO. Due to the 
slit curvature, spectra taken in B nodding positions are shifted 
in wavelength with respect to spectra taken in A nodding posi- 
tions. This was corrected for by mapping all ID spectra to the 
spectrum of the first A nod position. 

The individual CRIRES wavelength settings provide a con- 
siderable degree of spectral overlap and up to eight indi- 
vidual spectra were combined to one final spectrum at each 
wavelength. While merging the individual settings, small mis- 
matches in the wavelength solutions as well as imperfections 
in the individual spectra (detector cosmetics, ghost contamina- 
tion) were corrected. The final spectrum was normalized to a 
pseudo-continuum level of unity and fin ally shifted to match 
the McMath FTS spectrum of solar umbra (Wall ace et"aT]|1998l) . 
The error in the wavelength calibration should be smaller than 
0.75 km s" 1 . 

We find the SNR at the continuum level of most parts of the 
final spectrum to exceed 200. The high signal-to-noise ratio and 
high spectral resolution of the CRIRES data allow us to identify 
FeH lines with high accuracy. 
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Fig. 1. top: FeH Wing-Ford band, observed (black) and com- 
puted (red). Bottom: a magnification of the spectrum above. 



2.2. Theoretical FeH molecular data and line synthesis 

The theoretical data which we use to identify the FeH lines is 
taken from Dulick et al. (2003). They provide tables of quantum 
numbers and energiefl In particular, they provide the vibrational 
assignment of the upper and lower states v„ and v/, respectively, 
the projection of the total orbital angular momentum on to the in- 
ternuclear axis for the upper and lower state, Cl u and £2/, and the 
rotational quantum number// for the lower state. Furthermore 
the transition branch (P, Q, R), the parity, the wavenumber in 
crrT 1 , the lower state energy of the transition E/ and the Einstein 
A value are given. 

From this informati on it is possible to estimate \og l0 gf- 
values through (Bernath 2005) 



e G m e c . 

8 if = ^ — T1 A 8»' 
Lne L v L 



(1) 



with gi - 2Ji + 1 and g u - 2J U + 1 as the lower and upper 
statistical weights of the transition and v = c/A as the transition 
frequency. All quantities are in SI units. 

The van der Waal s broadening is determined following 
ISchweitzer et alJ ([1996) which is basically Unsold's hydrogenic 
approximation. For the ionization energy needed in this approx- 
imation we use an empirically determined value of 6 eV which 
was deduced from comparison with other diatomic molecules 



See http : //bernath . uwaterloo . ca/FeH 
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(IWendeetal.ll2009t) . Despite the fact that this is not the theo- 
retical value, which is slightly higher, we use this one since its 
influence on the van der Waals broadening is not significant. 

The molecular partition function for FeH QFeH, which 
is needed for the co ncentration of FeH is computed after 
ISauval & Tatum ( 1984 ) with m olecular data taken from Tables 9 
and 10 of Duli ck et al.l y003). We give in equation[2]a polyno- 
mial expression of a fit to the partition function which is valid 
between 1 000 K and 8000 K. 



QFeH = /_j aiT '> 



(2) 



where T is the temperature in K and 

( -4.9795007e + 02 ^ 
6.5460944e-01 
3.4171590e-04 
2.7602574e - 07 
1.0462656e- 11 ) 
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are the coefficients of the polynomial. For the creation of FeH 
i.e. the concentration, which is governed by the Saha-Boltzmann 
equation, we assume F e + H — > FeH. Wit h this partition func- 
tion and the data from Duli ck et aD (|2003), we can use a sim- 
ple description for the absorbance (described in section 13.41 in 
this paper in more detail) to compute FeH spectra with a sim- 
ple reversing layer model and separate the lines into vibrational 
bands in the observed wavelength region (see Fig. |2j). From this 
figure we can expect that there will be two sequences of vibra- 
tional transitions present. The sequence with Av = 0, which 
are the (0,0), (1,1), and (2,2) vibrational transitions, and the 
sequence with Av = 1, which are the (1,0), (2, 1), (3,2), and 
(4, 3) transitions. We note that this method does not take into ac- 
count (among many other things) the atmospheric structure or 
the chemical composition of the star. It gives only a rough esti- 
mate of the relative strengths of the FeH bands. 

The full synthetic line formation for the comparison and 
identification of the obs erved FeH lines i s done with the line 
formation code SYNTH3 (Kochukhov 2007). This code is able to 
compute large spectral regions with all FeH lines taken from our 
line list (see Fig.Q]). We use only lines with log 10 g/ > -7 be- 
cause we assume that other lines have no significant influence. 
By computing all the lines in a certain region at once, we account 
for blends, but we also compute all lines of the region individ- 
ually in order to measure their equivalen t widt h Wa using the 
LINF0R3D code (based on Basc heket alj|1966l). For the input 
model atmospheres we use MARCS (Gustafs son et alj |2008) with 
solar composition dGrevesse et alJl20O7l) . These models are well 
suited for these cool temperatures in low mass stars, since they 
make use of up-to-date atomic and molecular data and go down 
to effective temperatures of 2500 K. We use the plane-parallel, 
LTE models where the convection is treated in the mixing-length 
approximation. In the computation of these model atmospheres, 
the microturbulence parameter was set to zero. However, in the 
computations of the actual spectra, we assume micr oturbulence 
parameters according to the results of I Wende et alj (2009). We 
do not expect any significant influence since they are on the or- 
der of a few hundred m s _1 . We neglect the broadening from 
macro-turbulent motion, which would be hardly visible in the 
observed spectra. 

3. Methods 

We start with the investigation of the observed spectrum, for 
which we determine the position of the spectral lines and de- 
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Fig. 2. FeH vibrational bands separated by their vibrational 
quantum numbers. See color version for more details. 

cide whether a line feature is a blend or not. Then we measure 
the equivalent width W,\ of the lines with a Voigt fit procedure 
described below. 

We compare the line positions found in the CRIRES data 
to theoretical ones and identify them with FeH lines. In order to 
confirm an identification, we use statistical means: (i) the method 
of coincidence, and a cross-correlation technique producing co- 
incidence curves; and (ii) a method which takes the intensity into 
account. In this latter method we will compare theoretical line 
strength S (Honl-London factor) with observed Wa following 
Schadee ( 1964). We also use a description for the absorbance of 
spectral lines in order to correct theoretical line intensities given 
in terms of the Einstein A values. For this we will compare ob- 
served and computed spectra with each other and obtain a scal- 
ing factor for the Einstein A values. The final result is a corrected 
line list that reproduces the observed stell ar spectrum as well a s 
the line positions in sunspot spectra from Walla ce et alj ([1999). 
The line intensities are hard to confirm in the solar case, since 
many FeH lines are strongly split by magnetic fields. 

3.1. Voigt Fit 

In order to measure Wa in the observed spectra we use a 
'multi- Voigt fit' procedure (based on IDL curvefit function). The 
'multi- Voigt fit' is defined as 



I>, 



H(u<, ad 
max(H(ui, a,)) ' 



(4) 
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where A is the amplitude describing the depth of the line, 

A-A 

u - , 

<x 

with <x as Gaussian (or Doppler) width, 



An c cr' 



(5) 



(6) 



which will be called the Voigt constant throughout in this paper, 
y is the radiation damping constant, and H(u, a) is the Hjerting 
function (Gray 2008, and references therein). We successively 
fit the whole spectrum within bins of 7 A . The first and the last 
1 A are cut off after the fit to avoid boundary effects and we use 
only line profiles whose centers are inside the 5 A bin. For the 
Voigt profiles inside a bin, we assume a constant a-value for the 
Lorentz part and a fixed Gaussian width, but we allow changes 
in these two parameters from one wavelength bin to another be- 
cause the FeH lines tend to become narrower towards longer 
wavelengths. This is probably because we use a and <x only as fit 
parameters, and FeH lines starts to saturate at the band head, but 
become weaker towards longer wavelengths. Hence, the width 
of the line profile, which is a combination of a and cr decreases 
with decreasing saturation and so both parameters decrease. The 
wavelength dependence of the Doppler width affects the width 
of the lines as well, but it is negligible and goes in opposite di- 
rection (for example, cr - 0.1 A at A - 10000 A would change 
tocr = 0.105Aat/l= 10500A). 

With this fit procedure, we obtain the parameters (position, 
depth (amplitude), Voigt constant, and cr), needed for the indi- 
vidual Voigt line profiles to fit the observed spectrum. In Fig. [3] 
we show an example of how the fit (red) reduced the observed 
spectrum (black) into single line profiles (bottom panel). 

The convergence criterion is the minimization of the residual 
flux between observation and fit (O - C). The iteration is stopped 
if the maximal error of the fit is lower than three times the stan- 
dard deviation of the error or if the standard deviation of the error 
does not change significantly between two iterations. From these 
Voigt profiles, W\ can easily be computed by integrating over the 
single line profiles. The fit is also able to find and separate pos- 
sible blends. For this, we assume that the blended components 
differ in position by at least 0. 1 A . 

The measured Wa will be assigned to the associated theo- 
retical lines. This means if an observed line can be identified 
with exactly one theoretical line, Wa is fully assigned to this one 
theoretical line. However, in most cases theory predicts more 
than one line for an observed line position because of numerous 
overlapping vibrational bands in the same wavelength region and 
from many closely-spaced line pairs which only differ in parity. 
In these cases, Wa will be distributed to all predicted theoretical 
lines at this position. This is done using the ratio /?, of theoretical 



equivalent widths W\ obtained from the individual theoretical 
line profiles: 



iht't 
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W\ 



ZieB tke °W\ 



(7) 



with B as the set of blended components. There is of course 
the chance that FeH is not the only contributor to the observed 
blend. To avoid this uncertainty, we additionally use the line in- 
tensities in a method described below. 
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Fig. 3. Upper panel: Part of the observed spectrum of GJ1002 
with multiple Voigt fit. Lower panel: Single Voigt functions 
which were returned by the fit. 

3.2. Method of coincidence 

In order to confirm the identification of a molecular band, or a 
sub-band, we use the 'method o f coinc idence', which was first 
introduced by Russell & Bowenl (1 19291) . It gives the number of 
lines C in a spectral range that will be found by chance if one 
uses a set of randomly generated line positions and compare 
them with observed ones. Fundamental probability calculations 
lead to 



C = M[1 -exp(-2xw)], 



(8) 



where M is the number of lines in a particular region, x is the 
tolerated deviation in position, and w is the line density (average 
number of spectral lines in the investigated region). This means 
that C/M = ^random gives the probability of finding, with a ran- 
domly chosen line position, a random coincidence. If one iden- 
tifies N out of M lines in the observations, then the probability 
of finding an identification is N/M - /^identified- The ratio N/C 
describes how likely it is to find per identified line a line by 
random coincidence. Hence, the number of actually identified 
lines N should exceed the number of C of purely random coinci- 
dences. If this is the case, then one can assume that the lines are 
probably identified in the observed spectra. 

3.3. Theoretical Line Strength 

For the identification of molecular lines, it is also useful to take 
the intensities of the lines into account. In order to compare the 
theoretical lin e strength S w ith the observed equivalent width 
Wa, we follow lSchadeel (fl964). Forweak lines (mildly saturated) 
Wa is proportional to the wavelength Aq, the oscillator strength 
/, and Af;, the number of absorbers (Gray 2008): 



W A oc AlfNi, 
and 



81 



(9) 



(10) 



with Ni as the number of absorbers given by the Boltzmann dis- 
tribution 



Ni oc gie ( # E °\ 



(11) 
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where Eq is the lower state energy level in cm . In order to 
derive the final form, we use equation \W\ and QT| together with 
equation [9] and the expression for the Einstein A value 



^l',,l'/« 



(12) 



Here A VmVl is the Einstein A for a specific vibrational transition 
and constant for the vibrational band. The resulting equation is 



Wi 



he 



lo gi()7m = c - TT £ o l0 



SA$ 



kT 



(13) 



which is the relation from ISchaded (1 19641) with Eq replacing 
BJi(Ji + 1). Possible stimulated emission can be neglected for 
our analysis since there is only a small population in the excited 
state. If we plot log 10 Wa/S A^ against E , a straight line should 
be found, which then suggests the correct identification of the 
lines in the FeH molecular band. Equation [T3] accounts for dif- 



N(z) 



a and/? are free parameters and here chosen as 1.25 and -0.25, 
respectively. We plot the right hand side of equation [13] for con- 
stant and variable molecule number in Fig. [4] However, the sit- 
uation is much more complicated and we use synthetic line for- 
mation to investigate this behavior in the results section. 

Following Schadee (1964), we can now use the 
[log 10 Wa/S,Eq] diagram to classify the identified FeH 
lines into one of the following classes: 

1 . P - the line is present, and its Wa agrees well with the straight 
line of the diagram. 

2. Pb - the line is present, but its Wa is too large, i.e. the 
log 10 Wa/S value lies above the straight line of the diagram. 
This could imply that the line is blended by an other element 
(or that its computed line strength is too small). 

3. R - the line strength is presumably reduced by perturbations. 
That means, that the computed line strength is too large, and 
the data point lies below the line. 

4. Q - we identified the line, but we can't verify its identifica- 
tion, because we only investigate lines with AQ = in this 
plot. 



- 0.0010 - 



Eventually, it should be possible to derive the excitation tem- 
perature for the rotational transitions T in equation [13] from the 
slope of a linear fit i n the [ l ogm W a/S, Eq] diagram. However, 
one has to be careful. IWohll (Q/970) reports, and we experienced 
difficulties, that the obtained rotational temperature from this 
method relies crucially on the data points which are included 
in the linear fit and also on the degree of accuracy in measuring 
W A . 



1500 
J(J+1) 



Fig. 4. Relative number of absorbers with (dashed line), and 
without (solid line) constant number of FeH molecules over dif- 
ferent atmospheric layers. 

ferent branches as well as for different Q. This is because S on 
the left side of equation [T3l is computed for all branches (P, Q, 
R) and Q values. 

The S values that we will use in the analysis are computed 
for the intermediate Hund 's case and are dete rmined from the 
Einstein A values given bv lDulick et al.l d2003l) . 

We point out that we implicitly assumed in equation [TTI that 
the total number of FeH molecules N is constant for all lines. 
This assumption is only valid if we consider a small isother- 
mal atmospheric layer. However, since a spectral line forms over 
several layers, this assumption is not exactly valid. If we further- 
more consider a set of lines, with a wide range in gf and Eo 
values, then we expect that these lines form at different heights. 
Thus, the number of FeH molecules is not constant anymore: 
strong spectral lines are assumed to be formed in higher atmo- 
spheric layers, and weak lines in deeper layers. Hence, devia- 
tions from a straight line in the [log l0 Wa/S,Eu] diagram for 
lines with very small and very high J are expected. Assuming 
more FeH molecules in deeper layers, due to higher density, 
larger equivalent widths for these lines can be expected. For 
a qualitative description, we assume the total number of FeH 
molecules to be inversely proportional to the equivalent width, 
which reflects the heights of formation for weak and strong lines. 



N(z) = N ■ aW x . 



3.4. Line strength correction 

In some cases Wa and line depths of the observed FeH lines 
do not match the computed ones from the line list. In general, 
we observe that the differences increase towards longer wave- 
lengths and computed lines become stronger than observed ones. 
The strength of the lines are mainly determined by the lower 
state energies Eq and Einstein A values. We will correct only the 
Einstein A values, since a set of high resolution spectra in the z- 
band for different temperatures would be required to correct Eq. 
In order to corr ect the Einstein A value, we use the formula from 
iBernafhl d2005l) for the absorbance in a modified form in which 
the Einstein A values enter the expression 



*\{ 



(2J U + 1)A , __ 
8nv 2 q 



£o7«7l 



,-hv/kT 



)GNl, 



(15) 



where v is the frequency of a molecular in cm ~ ' , N is the number 
of molecular absorbers per cubic cm in the energy state (popu- 
lation density), q is the partition function, 1 is the length, e.g. for 
an atmospheric layer, and G is a line profile function, e.g. a Voigt 
function. If we compare the observed spectra with the computed 
one (produced with the SYNTH3 code) and assume that both at- 
mospheres have the same structure, then A x and E^ are the only 
variables that account for differences in the spectra (x stands for 
either the computed or observed spectra). We assume also E^ as 
constant and write 



(14) 



-Mr) 
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This shows that the correction for the Einstein A value is just a 
linear scaling with the ratio of the intensities 



-Mr) 

\ 'o Jobs 

-Mr) 

\ '0 ) comp 
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We introduced the scaling factor s for the r atio of both intensi - 
ties. If we use that A = A Vu>n S /(2J U + 1) dDulicketal.ll2003h . 
then 



t obs 



* obs 



(18) 



where S is the Honl-London factor. If the theoretical S values are 
wrong, we could perhaps detect it using equation [13] or perhaps 
in a possible dependence of the scaling factors on //. 

4. Results 

4.1. Atomic line identification and unidentified lines 

Since atomic lines are present in the wavelength region of the 
observed spectra of GJ1002, we use the VALD Q (Kupka et al. 
2000) database to include the available atomic line data in our 
calculations. We give a list of atomic lines which are present at 
Teg - 3100K in Table Q] We included lines with line depths 
deeper than 2% below the continuum in the computed spectra. 
We do not try to correct the atomic line positions or line strengths 
and take the data as provided by VALD. 

After we identify the atomic lines and FeH lines, there are 
still unidentified lines which seem to belong to neither FeH nor 
to a known atomic line. We give a list of these lines for which the 
line depth is deeper than 10% below the continuum (Table |2j. 
Our opinion is that most of these lines belong to FeH but we 
were not able to identify them with confidence. 

4.2. FeH line identification 

We use the results from the Voigt fit to assign the individual 
Voigt profiles from the observed lines to the individual theoreti- 
cal FeH lines. We did this by hand, and define a line as identified 
if the position of the observed line agrees within 0.1 A to the 
theoretical predicted position. If the offset is larger than 0.1 A 
we identify a line feature if the characteristic shape is similar 
in the observations and computations (e.g. the line at 999.55 nm 
or 999.8 nm in Fig. O. In Fig. [5] and also in the full FeH atlas, 
we labeled all identified lines with their quantum numbers, i.e. 
vibrational assignment (v«, v/), branch (P, Q, R), lower state ro- 
tational quantum number J, Q, and in the case of blends, with 
their contribution to the blend. We did this for the observed 
wavelength region and the complete plot is available in the on- 
line material of this paper. In Fig. [6]we show a histogram of the 
residuals between the computed line positions and the observed 
ones. There is no obvious systematic behavior of the scatter with 
wavelength (see inset in Fig.[6]l. The scatter follows a normal dis- 
tribution centered at AA = 0.02 A corresponding to 0.67 km s _1 . 
This mean value is beyond our spectral resolution and in gen- 
eral all residuals smaller than 0.75km s -1 (~ 0.025 A ) are not 
significant since they are also smaller than the accuracy of the 
wavelength calibration. From the investigation of the line posi- 
tions, we saw that the fraction of lines for which the residuals 
are smaller than a certain range are distributed as shown in the 
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Fig. 5. Observed spectrum of GJ 1002 (black) and computed one 
(red) labeled with quantum numbers. 

right inlay of Fig. [6] This plots shows that ~ 80 % of the line fea- 
tures could be identified by their positions which do not deviate 
by more than 0.1 A from the predictions. In the cases where the 
residuals are larger than 0.025 A the uncertainties in the molec- 
ular constants which were used to compute the FeH line list are 
responsible for these deviations. 
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Fig. 6. Histogram of the residuals between computed line po- 
sitions and observed ones. In the upper left corner, the resid- 
uals are plotted against wavelength. In the upper right corner, 
the fraction of lines with residuals lower than a certain value 
are plotted. The left dashed line represents the accuracy of the 
wavelength solution which is 0.025 A the right one the detec- 
tion boundary of 0.1 A . 

For our analysis we use only lines with Wa > 2 mA, which 
approximately describes a line with 2 % relative flux and a 
FWHM of 0.1 A. We ignore lines with smaller contributions be- 
cause their intensities are similar to the noise level. However, 
only 167 lines out of 1359 have equivalent widths less than 
2 mA. 

4.2.1. Vibrational bands 

From Fig. [2] we expect that the dominant vibrational bands are 
(0, 0) and (1, 1) from the Av = sequence, and (3, 2) and (4, 3) 
from the Av = 1 sequence. In Fig. Q we present a histogram 
with the number of identified lines for each vibrational band and 
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Table 1. Identified atomic lines 



Ton 


Position[A ] 


Ei 


l°g#/ 


Ion 


Position[A ] 


Ei 


logs/ 


Ion 


Position[A ] 


Ei 


logg/ 


'Crl' 


9903.6226 


2.987 


-2.131 


'Fel' 


10116.787 


2.759 


-3.705 


'Til' 


10399.651 


0.848 


-1.623 


'Til' 


9930.0728 


1.879 


-1.580 


'Til' 


10123.668 


2.175 


-1.722 


'Crl' 


10419.476 


3.013 


-1.806 


'Til' 


9944.1036 


2.160 


-1.821 


'Til' 


10148.300 


3.148 


-0.910 


'Fel' 


10425.885 


2.692 


-3.627 


'Til' 


9951.7317 


2.154 


-1.778 


'Fel' 


10148.342 


4.796 


-0.177 


'Til' 


10462.915 


2.256 


-2.054 


'Crl' 


9951.7997 


3.556 


-1.129 


'Fel' 


10157.947 


2.176 


-4.225 


'Fel' 


10472.522 


3.884 


-1.187 


'Til' 


9967.4679 


1.053 


-4.108 


'Fel' 


10170.256 


2.198 


-4.114 


'Crl' 


10489.119 


3.011 


-0.972 


'Til' 


10000.700 


1.873 


-1.840 


'Til' 


10173.274 


1.443 


-3.465 


'Til' 


10498.990 


0.836 


-1.739 


'Til' 


10005.828 


2.160 


-1.124 


'Fel' 


10197.900 


2.728 


-3.589 


'Crl' 


10512.887 


3.013 


-1.558 


'Til' 


10008.403 


1.067 


-3.626 


'Cal' 


10201.981 


4.680 


-0.369 


'Fel' 


10535.121 


3.929 


-1.482 


'Til' 


10014.490 


2.154 


-1.284 


'Cal' 


10205.743 


4.681 


-0.199 


'Crl' 


10552.983 


3.011 


-1.976 


'Scl' 


10027.787 


1.865 


-1.286 


'Cal' 


10205.803 


4.681 


-1.102 


'Til' 


10554.648 


1.887 


-2.607 


'Til' 


10037.242 


1.460 


-2.227 


'Cal' 


10211.457 


4.681 


-0.039 


'Fel' 


10580.037 


3.301 


-3.137 


'Til' 


10051.581 


1.443 


-2.205 


'Cal' 


10211.561 


4.681 


-1.102 


'Til' 


10587.533 


0.826 


-1.866 


'Scl' 


10060.261 


1.851 


-1.479 


'Fel' 


10221.204 


3.071 


-2.760 


'Til' 


10610.624 


0.848 


-2.761 


'Fel' 


10060.397 


5.033 


-1.231 


'Fel' 


10268.031 


2.223 


-4.533 


'Fel' 


10619.630 


3.267 


-3.128 


'Til' 


10060.485 


2.175 


-0.894 


'Cal' 


10291.397 


4.624 


-0.265 


'Crl' 


10650.558 


3.011 


-1.613 


'Til' 


10062.662 


1.430 


-2.351 


'Fel' 


10343.720 


2.198 


-3.574 


'Til' 


10664.544 


0.818 


-2.007 


'Fel' 


10067.804 


4.835 


-0.288 


'Cal' 


10346.655 


2.933 


-0.408 


'Crl' 


10670.437 


3.013 


-1.489 


'Til' 


10069.273 


2.160 


-1.750 


'Fel' 


10350.802 


5.393 


-0.548 


'Crl' 


10675.063 


3.013 


-1.374 


'Til' 


10077.885 


1.067 


-4.065 


'Fel' 


10381.844 


2.223 


-4.145 


'Til' 


10679.972 


0.836 


-2.592 


'Crl' 


10083.115 


3.556 


-1.307 


'Til' 


10393.591 


1.502 


-2.595 


'Fel' 


10728.124 


3.640 


-2.763 


'Fel' 


10084.158 


2.424 


-4.544 


'Crl' 


10394.793 


3.010 


-2.006 


'Til' 


10729.329 


0.813 


-2.156 


'Crl' 


10114.770 


3.013 


-2.073 


'Fel' 


10398.645 


2.176 


-3.390 


'Cal' 


10729.654 


4.430 


-1.841 



Table 2. List of unidentified lines deeper then 0.9. 



Ivac [A ] 


Depth 


l mc [A ] 


Depth 


Kac [A ] 


Depth 


hac [A ] 


Depth 


hac [A ] 


Depth 


hac [A ] 


Depth 


/W [A ] 


Depth 


9904.46 


0.88 


10084.8 


0.90 


10235.9 


0.84 


10344.9 


0.86 


10424.9 


0.79 


10547.5 


0.80 


10644.2 


0.84 


9909.42 


0.87 


10095.7 


0.86 


10238.5 


0.73 


10346.0 


0.76 


10427.8 


0.89 


10548.2 


0.85 


10645.7 


0.84 


9927.74 


0.79 


10098.8 


0.89 


10242.9 


0.89 


10347.3 


0.82 


10438.8 


0.87 


10548.7 


0.89 


10655.3 


0.80 


9930.70 


0.81 


10107.3 


0.85 


10246.4 


0.81 


10347.5 


0.87 


10440.0 


0.85 


10550.8 


0.88 


10660.2 


0.79 


9931.26 


0.87 


10107.6 


0.86 


10247.3 


0.74 


10347.7 


0.87 


10442.0 


0.80 


10555.8 


0.90 


10660.7 


0.87 


9932.82 


0.88 


10108.7 


0.89 


10247.3 


0.74 


10351.3 


0.89 


10443.2 


0.84 


10559.0 


0.89 


10665.1 


0.84 


9932.94 


0.88 


10116.9 


0.65 


10248.4 


0.84 


10354.4 


0.79 


10443.6 


0.89 


10563.5 


0.82 


10665.6 


0.71 


9938.22 


0.88 


10140.2 


0.89 


10255.4 


0.84 


10357.5 


0.89 


10444.6 


0.90 


10563.8 


0.67 


10666.9 


0.87 


9978.82 


0.89 


10141.5 


0.73 


10260.6 


0.88 


10358.5 


0.74 


10444.8 


0.81 


10565.8 


0.82 


10669.6 


0.88 


9983.94 


0.90 


10144.5 


0.66 


10261.0 


0.88 


10361.5 


0.86 


10444.9 


0.90 


10567.1 


0.83 


10671.3 


0.90 


9984.94 


0.90 


10164.4 


0.89 


10266.4 


0.81 


10363.3 


0.89 


10446.7 


0.88 


10568.8 


0.90 


10671.6 


0.71 


9985.22 


0.90 


10167.5 


0.84 


10267.3 


0.79 


10367.9 


0.87 


10449.3 


0.88 


10571.1 


0.89 


10673.3 


0.88 


9993.22 


0.72 


10171.2 


0.89 


10274.1 


0.74 


10369.1 


0.88 


10463.2 


0.81 


10577.5 


0.90 


10674.7 


0.83 


9993.82 


0.77 


10174.1 


0.89 


10280.7 


0.87 


10378.0 


0.83 


10463.9 


0.84 


10580.6 


0.85 


10675.5 


0.81 


9994.14 


0.87 


10177.2 


0.90 


10281.0 


0.77 


10381.0 


0.88 


10464.3 


0.77 


10582.8 


0.80 


10681.4 


0.90 


9997.30 


0.88 


10187.0 


0.81 


10286.1 


0.89 


10382.0 


0.90 


10464.9 


0.90 


10586.5 


0.87 


10697.2 


0.89 


9999.94 


0.77 


10188.6 


0.67 


10287.4 


0.82 


10382.7 


0.78 


10477.6 


0.75 


10586.9 


0.75 


10697.3 


0.88 


10005.3 


0.89 


10193.7 


0.83 


10297.2 


0.80 


10383.8 


0.85 


10479.6 


0.86 


10588.5 


0.90 


10711.9 


0.81 


10015.0 


0.81 


10197.2 


0.87 


10298.0 


0.88 


10383.9 


0.85 


10490.9 


0.83 


10589.6 


0.89 


10712.3 


0.86 


10021.2 


0.88 


10200.1 


0.89 


10303.5 


0.88 


10385.7 


0.89 


10493.9 


0.87 


10593.8 


0.84 


10713.0 


0.74 


10023.7 


0.87 


10202.9 


0.90 


10309.6 


0.86 


10387.8 


0.86 


10499.5 


0.83 


10594.3 


0.76 


10716.6 


0.85 


10027.7 


0.89 


10203.0 


0.90 


10310.6 


0.84 


10392.2 


0.87 


10499.6 


0.83 


10600.7 


0.87 


10724.6 


0.77 


10029.3 


0.86 


10209.1 


0.86 


10312.2 


0.75 


10396.0 


0.73 


10500.2 


0.75 


10601.8 


0.89 


10725.3 


0.81 


10031.8 


0.70 


10211.1 


0.70 


10312.5 


0.89 


10398.4 


0.89 


10500.9 


0.87 


10602.9 


0.88 


10726.4 


0.86 


10035.9 


0.83 


10214.5 


0.79 


10314.7 


0.84 


10399.2 


0.87 


10514.7 


0.84 


10604.0 


0.89 


10734.7 


0.85 


10041.9 


0.90 


10215.8 


0.77 


10319.3 


0.90 


10400.2 


0.82 


10515.1 


0.88 


10605.3 


0.90 


10739.2 


0.90 


10048.2 


0.82 


10217.8 


0.89 


10319.5 


0.85 


10400.5 


0.88 


10515.5 


0.79 


10605.5 


0.88 


10740.5 


0.82 


10051.1 


0.86 


10219.5 


0.86 


10320.6 


0.86 


10402.8 


0.83 


10516.6 


0.89 


10609.5 


0.87 


10744.1 


0.86 


10066.5 


0.88 


10219.9 


0.87 


10322.2 


0.88 


10406.5 


0.89 


10519.3 


0.88 


10609.6 


0.88 


10745.3 


0.84 


10068.7 


0.79 


10220.0 


0.88 


10322.8 


0.75 


10406.8 


0.81 


10519.7 


0.89 


10612.7 


0.77 


10748.7 


0.90 


10073.9 


0.58 


10223.9 


0.85 


10327.8 


0.89 


10408.4 


0.87 


10532.9 


0.86 


10613.1 


0.82 


10749.4 


0.88 


10074.5 


0.89 


10225.8 


0.77 


10327.9 


0.87 


10409.3 


0.75 


10535.5 


0.89 


10614.0 


0.81 


10752.1 


0.90 


10075.9 


0.90 


10226.7 


0.76 


10328.3 


0.78 


10410.1 


0.81 


10536.2 


0.86 


10626.7 


0.82 


10753.7 


0.90 


10078.5 


0.86 


10227.7 


0.88 


10332.7 


0.88 


10415.7 


0.89 


10539.5 


0.83 


10629.9 


0.88 


10754.7 


0.79 


10080.6 


0.75 


10229.4 


0.89 


10338.9 


0.64 


10417.8 


0.88 


10543.1 


0.87 


10634.8 


0.89 


10755.8 


0.90 


10083.8 


0.72 


10230.5 


0.73 


10340.0 


0.89 


10422.5 


0.90 


10546.1 


0.86 


10638.3 


0.80 


10759.4 


0.89 
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Table 3. Results from the coincidence method for lines with 

W A >2mA. 



Band Wavelength [A ] C w[#/A] 



N(M) 



(0,0) 


9900- 


10762 


360 


3.73 


489 (658) 


(1,0) 


9901- 


10759 


52 


3.74 


25 (99) 


(1,1) 


9941- 


10759 


208 


3.71 


251 (388) 


(2,1) 


9931- 


10706 


68 


3.76 


45 (131) 


(2,2) 


10480 


- 10764 


65 


3.81 


26 (123) 


(3,2) 


9910- 


10764 


79 


3.78 


52 (150) 


(4,3) 


9905- 


10673 


119 


3.75 


68 (226) 



shows that the expected vibrational bands are present. We also 
show the number of possible lines from theory with W,\ > 2 mA 
in the observed wavelength region. This number is based on line 
by line computation of Wa- In this histogram the bars in the fore- 
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Fig. 7. Histogram of predicted and identified FeH lines in the 
GJ1002 spectrum. 



ground take into account only the lines identified by their posi- 
tions. 



4.2.2. Coincidence-by-Chance Method 

For the lines with residuals smaller than 0.1 A we compute the 
coincidence-by-chance factor C from equation[8] In Table[3]the 
values of C, w, N, and M are given for all identified bands. For 
the (0,0) and (1,1) transition, the number of identified lines N 
exceeds the number of coincidences by chance C, which is a 
clear indication that these bands are present in the investigated 
region. For the other bands the situation is not so clear and we 
will confirm them with further investigations. 



4.2.3. Cross-Correlation Method 

We use cross-correlation techniques to investigate the agreement 
between the theoretical line list and the observed spectra. As 
a reference and a test, we cross-correlate a computed spectrum 
from this theoretical line list, which is broadened by an instru- 
mental profile with a resolving power of 70 000 , with the origi- 
nal line list (e.g. lFawzvll 995: Faw zv et aD[l99 8 i ). To be specific, 
we vary the theoretical positions with steps of A<x = 0.0125 A 
in a range of 0.375 A and measure the relative intensity, at the 
different positions, weighted with Wj of the line and integrate 
over all lines. We then normalize the results with the number of 



lines in the line list. We do this for all lines in the vibrational 
bands which are found in the M dwarf spectra (see Fig. [8). If 
a vibrational band is present, a peak around zero appears above 
the noise produced by random coincidences with other lines. 

We produce three different curves. For the first one, we use 
all possible lines from the original line list for comparison with 
the observed stellar spectrum (solid line in Fig. [HJ. For the sec- 
ond curve we compute a reference curve by cross correlating 
the original line list with a synthetic spectrum computed from it 
(dashed line in Fig. |H). This case produces the maximum pos- 
sible correlation. A third curve is produced using the corrected 
line list containing only identified FeH lines and we compare it 
with the observed stellar spectrum (dotted line in Fig.[HJ. 

As expected, for the synthetic spectrum with the theoreti- 
cal line list (dashed line), all bands show clear peaks above the 
noise. For the observed spectrum (solid line) and original line 
list, the (0, 0) and (1,1) vibrational bands show clear peaks above 
the noise, which is in agreement with the coincidences by chance 
values in Table[3] After improving the FeH line list (dotted line), 
all bands show peaks above the noise similar to the peaks ob- 
tained in the reference case from cross -correlating the theoreti- 
cal line list with the computed spectra. The original theoretical 
line positions for the (1,0), (2, 1), (2, 2), (3, 2), and (4, 3) bands 
are not accurate enough to show significant peaks in the cross 
correlation. To confirm our identifications, we also use the line 
strength in the next section. 




synthetic model 

obs. before line correctio 

obs. after line correction 



-0.2 0.0 



Fig. 8. Cross-correlation curves for different vibrational bands. 
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Fig. 10. Logarithm of computed equivalent width and theoreti- Fig. 11. log 10 r mra „ and Wa are plotted as a function of lower 
cal line strength against lower level energies for the (0, 0) band, state energy for computed FeH lines from the (0, 0) band for the 
In each plot a linear fit to the data is shown. Different colors of three branches, 
the data points belong to different Q. values. Dots represent lines 
which are identified within 0.1 A and diamonds lines which dif- 
fer more than 0.1 A from the theoretical position. 



4.2.4. Line Intensity Method 



We use the method described in Secf3]to test if the tentatively 
identified observed FeH lines can be identified with the theo- 
retical ones. We plot log 10 Wa./SAz against Eq (see equation [TBI 
for the identified vibrational bands, branches, and Q., using an 
appropriate estimated error (see Fig. |9). 

In almost all cases a linear correlation is visible. In some 
cases where lines with small J are present, we see a deviation 
from the straight line. This deviation is expected due to differ- 
ent heights of formation for the FeH lines. We calculated model 
plots with synthetic lines, which reproduces this behavior in de- 
tail (Fig. [Toll. The computations are very similar to the observa- 
tions, which supports the assumption of different heights of for- 
mation. To test this assumption, we computed the contribution 
function for the FeH lines (Magain 1986) and use them to de- 
termine a weighted mean of the continuum optical depth, T mean . 
The latter describes a representative atmospheric line formation 
depth. In Fig. [TTJ we plotted W,\ and T mean against Eq for each 
branch from the (0, 0) band. These plots show clearly that for 
the P- and R-lines with low and very high Eq originate in deeper 
layers than the lines with medium Eq. We also see that the lines 
with low and high Eq are the ones with small Wa because the 
Honl-London factors for P- and R-lines are proportional to // 
and hence the lines become stronger for large //, but they also 
decrease for large J\ due to their increasing E$&. Hence, there is a 
maximum in Wa for medium J/. In the case of Q-branch lines, the 
situation is different. The Honl-London factors decrease with in- 
creasing J i and hence the line strength mono tonic ally decreases 
with 7/. In this case, the lines with low Eq are the lines with 
the largest Wa and hence these lines are formed in higher lay- 
ers and the Q-branch shows the opposite behavior to the P- and 
R-branches in the log 10 -^4- plot. 

We can conclude that in deeper layers, where weak lines are 
formed, the equivalent width is enhanced by the larger num- 
ber of FeH molecules that contribute to the absorption (see 
Fig. fT2l) . The molecule number increases towards deeper lay- 
ers due to the higher overall density even though the concentra- 
tion of FeH molecules relative to Fe and H decreases. Towards 




Fig. 12. The number of FeH molecules as a function of optical 
depth in the region where the FeH absorption lines are formed. 



higher temperatures in deeper layers or in hotter stars, the num- 
ber of FeH molecules would decrease again due to the ionisation 
of Fe which is then no longer available for the formation of the 
molecule. In our case, this behaviour results in a deviation from 
a straight line in Fig. [9] 

With this knowledge, we can use the line strength to con- 
firm the identification of FeH . In Fig. |9l clearest correlations 
are found for the (0,0) and (1,1) bands which have also the 
largest number of lines, but also for the P-branch of the (1,0) 
and (2, 1) transitions and the R-branch of the (2,2) transition. 
The (3, 2) and (4, 3) transitions show also linear correlations but 
with a larger scatter. This result is a strong indication that all the 
identified bands are present in the observed spectra. Although 
the Q-branches of the (0, 0) and (1,1) bands show no clear linear 
dependence for large Eq, we expect them to be present and cor- 
rectly identified, since they show the expected downward trend 
for small energies. The scatter towards higher energies is much 
larger than in the other branches and there are only a few mea- 
surements. The slope seems to be positive, but we could expect 
it to be negative if we would have more data points. 

We fit the data points linearly for lines with // > 7 (to avoid 
the region where the influence of variable FeH number is too 
strong) and use the difference of the data points from the fit as a 
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Fig. 9. Logarithm of observed equivalent width and theoretical line strength against lower level energies. In each plot, a linear fit to 
the data is shown. Different colors of the data points belong to different Q. values. Dots represent lines which are identified within 
0. 1 A and diamonds lines which differ more than 0. 1 A from the theoretical position. The dotted lines represent the three <x scatter 
around the fit 



S. Wende et al.: CRIRES spectroscopy and empirical line-by-line identification of FeH molecular absorption in an M dwarf 
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0.0384 
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P 



measure for the confidence of an identified line. A line is classi- 
fied as "P" if the difference between log, -^ and the linear fit 

is smaller than three times the standard deviation <x of the scat- 
ter around the fit. But we assume that lines with / < 7 are also 
present, since they show exactly the expected behavior. It is clas- 
sified as "Pb" if the deviation is greater than cr and the data point 
is above the linear fit and as "R" if it is below. It is classified as 
"Q" if the line could not be investigated because we only investi- 
gated lines with AO = 0. We give a list of all identified lines with 
quantum numbers and corrected wavelength in Table [4] This list 
is explained in more detail in the appendix. 



4.3. Corrections to the line strengths 

If we want to use equation [T7] to correct for the differences in 
line depth, we have to match the stellar parameters as closely 
as possible. These parameters are basically effective tempera- 
ture, surface gravity and chemical composition, as well as van 
der Waals broadening constants, whose influence become sig- 
nificant at these low temperatures. Van der Waals broadening is 
computed with Unsold's hydrogenic approximation, and an en- 
hancement factor is used model the line wings correctly. Since 
for FeH no enhancement factor is reported, we need to determine 
one. ForM-type stars, an assumed surface gravity of log g = 5.0 
is standard and the chemical composition is usually assumed as 
solar. To match the strong Ti lines in the 10 300-10 700 A re- 
gion as well as the Fe H lines, we have to in crease the iron abun- 
dance from 7.41 to the lGrevesse & And ers ( 1989) value of 7.63. 
We use this scaling as a parameter and do not claim this to be 
the actual iron abundance of GJ1002. The free parameters for 
GJ1002 are now r e ff, van der Waals broadening enhancement 
factor (which we call from now on p v dw), and the instrumental 
resolving power, which we use as a fitting parameter to account 
for possible additional rotational broadening. These three param- 
eters are strongly correlated. We create^- 2 maps to determine the 
most likely combination that matches the observed spectra best. 



4.3.1. x maps 

For the comparison between observation and computation we 
chose a region where the original line list fits best, the lines are 
strongest and hence the influence of van der Waals broadening 
is largest. We selected the first 100 A from the (0,0) band head 
at 9900 A . For the computations we use our new line list with 
corrected positions and include also the identified atomic lines. 
To create the x 2 maps for the three combinations of parameters, 
we search for the minimum for each parameter (light cross in 
left plots in Fig. [TBI and use this value to construct the x 1 map 
for the other two parameters. The x 1 maps (right hand side in 



Fig.fTJb yield a consistent picture of the parameter combinations 
for the spectra of G J 1 002 . 





3000 3050 



(150 5200 



3 ^H 

rflfl 




Fig. 13. Left: Minima in x 2 plots for the resolving power, T e g , 
and fivdw (from top to bottom). Right:^- 2 maps for the three pa- 
rameter combinations, j3 v dw - T e $, fi v dw - resolving power, and 
resolving power - T e g (from top to bottom). The x 2 maps are 
computed for the minimum value of the leftover parameter on 
the left hand side. 

The most likely parameters for effective temperature, resolv- 
ing power, and/3 v dw are 3100 K, 70000 and 1.75, respectively. It 
is obvious, however, that the x 2 curves show regions with broad 
minima which would allow for variations in the derived param- 
eters. 

Since the observed spectra are supposed to have a resolving 
power of 100 000, the difference from our determined resolving 
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power stems probably from rotation, which is measured to be 
lower then 3 km s , but not necessarily zero. The difference in 
resolving power results in v sin i a 1.3 km s at a wavelength of 
10 000 A. 

An independent constraint for the instrumental resolving 
power and effective temperature is given by the Ti lines, which 
are strong and distributed over a wide wavelength range in the 
spectra. The computation of these lines with the parameters 
and broadening constants given by VALD fits the observations 
within 5 % for the line depth. This gives us confidence that the 
decreasing line depth of the FeH lines with increasing wave- 
length is a real feature and not due to normalisation effects. 

4.3.2. Einstein A values 

In order to correct for the Einstein A values we use computed 
spectra which have already corrected line positions. We itera- 
tively adjust the Einstein A values because for saturated lines 
the first scaling is not sufficient. The scaling factors for each line 
are listed in the Table [4] and an example for the corrected com- 
puted spectra is shown in Fig. [14] In order to estimate an error, 
we assume an accuracy of 1 % for the observed line depth, which 
results in an error of ~ 3% for the scaling factors. If we fur- 
thermore assume that the line depth is modified by an unknown 
blended feature by, e.g., 5% an error of ~ 16% follows. The ac- 
curacy of the Einstein A scaling is also influenced by the van der 
Waals enhancement factor f3 v dw- A change of +1 gives a mean 
difference of ~ 5% in the scaling constants, but can be up to 30% 
for individual lines, due to the logarithmic ratio of the intensities 
(see e.g. equation[T7l>. 

We obtain a good fit to the data with the scaled Einstein A 
values even for some lines calculated to be very weak. However 
this results in some cases in unrealistic log gf values, and we 
assume that these weak lines are blended with unknown compo- 
nents. The scaling of line blends is a difficult problem because 
it results in equal scaling factors for lines with completely dif- 
ferent quantum numbers. To avoid this problem one could deter- 
mine scaling factors for each branch, but we chose the simpler 
scheme of scaling each line. We use equation [T7] and plot s for 
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Fig. 14. Observed spectrum of GJ1002 (black) and computed 
one before A correction (red) and after A correction (green) la- 
beled with quantum numbers (both with corrected positions). 



the Einstein A values against 7; to search for a possible rotational 
dependence (Fig. [T3T >. We plot only scaling factors for the (0, 0) 




Fig. 15. Scaling factor for the Einstein A values against rota- 
tional quantum number for the (0, 0) and (1,1) bands. The plot 
is truncated at J/ = 35 for better visibility. 



and (1,1) bands, since these are the bands with the largest num- 
ber of lines and we use only lines which contribute more than 
99 % to a blended line feature in order to avoid contributions 
from incorrectly scaled lines. 

We find two groups of scaling factors. One group describes 
a strong enhancement of the line depths ('positive' scaling fac- 
tors) and the other only a small enhancement for low J\ and a 
reduction towards lines with large 7/ ('negative' scaling factors). 
We divide these two groups by a dashed line in Fig. [15] For the 
AQ = transitions of the (0,0) and (1,1) bands we include lin- 
ear fits to the data in Fig.[l5]to indicate the slope. 

The group of positive scaling factors is strongly dominated 
by AQ+ 1 transitions (diamonds and rectangles), while the group 
of negative scaling factors consists of only AQ = transitions 
(circles). The latter scaling factors for (0,0) and (1,1) bands 
(black and red circles, open for P-branches and filled for Re- 
branches) show an almost linear behaviour with J\ and become 
stronger towards larger J\. The positive scaling factors show also 
two linear groups, which originate from the (1, 1) band (red rect- 
angles) and from the (0, 0) band (black diamonds and rectan- 
gles). These scaling factors describe a strong enhancement of 
the lines for small 7/ and become smaller towards large J\. All 
groups of scaling factors have a similar negative slope, and are 
only shifted by a constant factor to higher or lower scaling fac- 
tors. 

The distribution of other bands and lines which contribute 
less than 99 % to a blended feature gives only a larger scatter 
to the data points, but does not change the basic trend of the 
scaling factors. We conclude that the J\ dependence in the scal- 
ing factors likely indicates shortcomings in the calculated Honl- 
London factors. In particular, satellite branches with AQ + 1 are 
much stronger than expected for Hund's case (a). In other words, 
the two 4 A electronic states are heavily mixed with other elec- 
tronic states and the simple Hund's case (a) behaviour antici- 
pated for isolated electronic states with relatively large spin-orbit 
splittings is not found. 



4.4. Rotational temperatures 

A rotational temperature T ro , can be obtained from the slope m 
of the linear fit in the [log 10 Wa/S, Eq] diagram described in the 



S. Wende et al.: CRIRES spectroscopy and empirical line-by-line identification of FeH molecular absorption in an M dwarf 



13 




Fig. 16. Rotational temperatures derived from the slope of the 
linear fits in Fig. [9] The error bars indicate the one sigma level. 
The dotted lines give the expected upper and lower effective tem- 
peratures for a 5.5 M-dwarf. The grey solid line is the weighted 
mean of the rotational temperatures with its one sigma error 
(dashed-line). 



above section. T rot can be calculated from 
he 



(19) 



For the fit to the data, we use only lines with 7/ > 7 in order 
to avoid the significant influence from varying absorber number, 
and neglect lines with line depth greater than 0.5 to avoid satu- 
ration effects. 

Due to the large errors (derived from the uncertainty of the 
slope of the linear fit in Hig- [9j, we consider only rotational tem- 
peratures with moderately small one sigma errors (see Fig.[T6|). 

We find systematically lower temperatures for the P- 
branches in comparison to the R-branches which is consistent 
with different heights of formation for most of the lines in a 
branch as seen in Fig.fTTI 

If we compute the weighted mean of the rotational tempera- 
tures (grey solid line in Fig. [Toll with the lcr-error (grey dashed- 
line in Fig. [Toll, we obtain T ro , « 3200 ± 100 K . This is in the 
middle of the expected temperature range for this spectral type 
(3000K-3300K, dashed-dotted lines in Fig.[[6]l and is close to 
our estimated value of ~ 3100K. The main contribution to T rot 
stems from the P- and R-branches of the (0, 0) transition due to 
their small uncertainties. The weighted mean of the P- and R- 

branches are Y rot * 2600 ± 150 K and T*, « 3750 + 150 K, 
respectively. 

We want to point out that even for a single branch in a band, 
the obtained temperature is an average over the individual exci- 
tation temperatures for each line. The resulting temperature de- 
pends crucially on the selection of the lines that are used. If one 
uses lines with similar equivalent width and lower level energy, 
then one could obtain temperatures for certain regions in the at- 
mosphere. For this method, however, a large number of lines is 
required, because otherwise the uncertainties become too large. 

We conclude that in order to use the method of rotational 
temperatures, a large number of well measured lines are required 
to minimize the error in the slope. Finally, the rotational temper- 
ature can only be expected to match the effective temperature if 
the lines form in a region around optical depth unity. 



5. Summary and conclusion 

We investigated the z-band region of the M5.5 dwarf GJ1002 
with high-resolution CRIRES spectra. This is the region where 
the (0, 0) vibrational band of FeH is present ('Wing-Ford band'). 
We were able to identify the (0,0), (1,0), (1, 1), (2, 1), (2,2), 
(3,2), and (4,3) bands. For confirmation of the band assign- 
ment, we used the method of coincidence, cross-correlation tech- 
niques, and the line intensities. 

For the identified lines, we applied empirical corrections to 
the theoretical line positions. The deviations between observed 
and computed positions are Gaussian distributed around zero. 
In the case of small deviations (< 0.025 A ), this could be due 
to uncertainties in the wavelength calibration of the observed 
spectra and in other cases due to uncertainties in the molecular 
constants which were used to generate the line list. Note that 
these empirical wavelength adjustments reproduce the observed 
stellar spectrum, but the quantum number assignments are by no 
means assured. 

The method of coincidence confirms the presence of the 
(0, 0) and (1,1) bands, but for the other bands we needed cross- 
correlation techniques to show their presence. Again the (0, 0) 
and (1,1) bands show clear peaks in the cross-correlation func- 
tion even with an uncorrected line list, but for the other bands 
it was necessary to use corrections to the line list in order to 
confirm them. We also used line intensity information by plot- 
ting the ratio W,\/S A A Q against Eq which shows a linear behaviour 
if the band is present in the spectra. With this method we could 
confirm the presence of all the other bands although not uniquely 
for the Q-branches. 

From the slope of the line in the [Wa/S Aq, Eq] plots, it was 
possible to derive excitation temperatures for rotational transi- 
tions, which could be identified with the effective temperature of 
the star if the lines are formed in the photosphere. We showed 
that this method is very uncertain since the error in the slope 
is high. The derived temperatures for the individual vibrational 
transitions range from ~ 2500 K to ~ 5500 K for GJ1002, but 
the weighted mean T ro , ss 3200 + 100 K is very close to the ex- 
pected temperature of an M 5.5 dwarf. However, the large error 
bars and differences between P- and R-branch temperatures sug- 
gests that this agreement may be more of a coincidence than a 
physical result. 

Finally we corrected the line strength in the FeH line list 
by scaling the Einstein A values, as some of these lines show 
large discrepancies compared to the observations. For this pur- 
pose it was necessary to derive the instrumental broadening 
(which included the rotational broadening), effective tempera- 
ture, and an enhancement factor for the van der Waals broaden- 
ing constants. The instrumental resolving power was derived to 
R - 70 000, which is equivalent to a rotational broadening with 
vsin(i) ~ 1.3 km s -1 . We also derived an effective temperature 
of r e ff = 3100K, and a van der Waals enhancement constant of 
1.75. The scaling factors of the Einstein A values show an al- 
most linear dependence on //, which indicates that there is likely 
a problem in the calculation of the Honl-London factors. 

With the improved identification of FeH lines, it is now pos- 
sible to characterize the FeH lines in the z-band region (e.g. 
magnetically sensitive and insensitive lines, temperature sensi- 
tivities of individual lines). Our improved line list will aid in 
the identification and simulation of FeH lines in spectra of cool 
stars. 
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Appendix 

Explanation of the FeH Table 

The Table|4]contains the observed wavelength A„b s which are ob- 
tained with the Voigt fit and the theoretical wavelength A,i, eo from 
the list of Dulick et al. (2003). The wavelengths are in vacuum. 
We give also the quantum numbers of the lines and the Einstein 
A values with their scaling factors sa- The lower level energy 
Ei is given in eV. The difference in position AA = A,/ wo - A \, s 
is also printed in the table. If the line is a blended line, then its 
contribution to the blend is given as the fraction normalized to 
one. If the line is not blended its blend value is one. We give 
then the classification of the line, as defined in section 3.3. We 
add a comment if the line is blended by an atomic feature, or if 
the classification of the line did not agree with the scaling factor 
of the Einstein A values. 
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Explanation of the FeH Atlas 

We plotted the whole spectrum in bins of 1 nm from 990 nm to 
1076.6 nm. Shown are the observed spectrum of GJ 1002 (black), 
the computed spectrum with corrected positions (red), the com- 
puted spectrum with corrected positions and scaled Einstein A 
values (green). We also labeled all lines with W,\ > 2 mA with 
quantum number for the vibrational transition, the branch, the 
lower J, the upper and lower Q, and in the last position, their 
blend fraction. The blend fraction is unity if a line is not blended. 
We also labeled the position of atomic lines with the element 
name below the spectrum. 
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